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Abstract. The instrumental line background of the Compton telescope COMPTEL onboard the Compton Gamma- 
Ray Observatory is due to the activation and/or decay of many isotopes. The major components of this background 
can be attributed to eight individual isotopes, namely 2 D, 22 Na, 24 Na, 28 Al, 40 K, 52 Mn, 57 Ni, and 208 T1. The 
identification of instrumental lines with specific isotopes is based on the line energies as well as on the variation of 
the event rate with time, cosmic-ray intensity, and deposited radiation dose during passages through the South- 
Atlantic Anomaly. The characteristic variation of the event rate due to a specific isotope depends on its life-time, 
orbital parameters such as the altitude of the satellite above Earth, and the solar cycle. 

A detailed understanding of the background contributions from instrumental lines is crucial at MeV energies 
for measuring the cosmic diffuse gamma-ray background and for observing 7-ray line emission in the interstellar 
medium or from supernovae and their remnants. Procedures to determine the event rate from each background 
isotope are described, and their average activity in spacecraft materials over the first seven years of the mission 
is estimated. 

Key words. Methods: data analysis - Line: identification 



1. Introduction 

Gamma-ray experiments in low-Earth orbit, such as the 
Compton telescope COMPTEL onboard the Compton 
Gamma- Ray Observatory (CGRO), operate in an in- 
tense and variable radiation environment. The main con- 
stituents of the ambient radiation fields are primary 
cosmic-ray particles, geomagnetically trapped radiation- 
belt particles, as well as albedo neutrons and 7-ray pho- 
tons. The different particle species interact with the space- 
craft and detector materials, resulting in the emission of 
instrumental background photons (for a review, see e.g. 
Dean et al. |l99l| ). COMPTEL data, dominated by instru- 
mental background, have a typical signal-to-noise ratio of 
a few percent. Hence, a qualitative and quantitative under- 
standing of the instrumental background is crucial for con- 
ducting astrophysical measurements, in particular of the 
cosmic diffuse gamma-ray background (hereafter CDG), 
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and of the 7-ray line emission in the interstellar medium 
or from supernovae and their remnants. 

The instrumental background experienced by 
COMPTEL is subdivided into two major compo- 
nents according to their signature in energy space: first, 
a continuum background discussed by Ryan et al. ( 1997 ); 
second, the instrumental line background, the focus of 
this paper. The latter arises from a number of different 
radioactive isotopes generated in the instrument material. 
This primarily occurs from activation by trapped protons 
during passages through the South-Atlantic Anomaly 
(SAA), from neutron absorption, and from primordial 
radioactivity^ An earlier report on activation in the 
COMPTEL telescope was given by Morris et al. ( |l997a| ). 

The discussion is structured as follows. After a brief 
description of the COMPTEL instrument in Sect. ||, gen- 
eral characteristics of the instrumental (line) background 
are summarized in Sect. [}| In Sect. |> identifications of 
specific isotopes are discussed. In Sect. |, the variations of 
the activity of individual isotopes are described. In Sect. [}| 
a comparison of instrumental line backgrounds in different 
low-energy 7-ray experiments is given. Also, average val- 



In the following, "primordial radioactivity" denotes unsta- 
ble isotopes that already existed when the Earth was formed. 
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ues for the activity of spacecraft materials are presented. 
The results of this work are summarized and discussed in 
Sect. |7[ Finally, appendices give the event selections used 
in these line studies, and provide detailed descriptions of 
the procedures employed for determining the background 
contributions of individual isotopes in the CDG analysis, 
and also - with slight modifications - in the analysis of 
the galactic 1.8 MeV line emission from 26 Al. 

2. Instrument description 

COMPTEL is the first double-scattering Compton tele- 
scope designed for 7-ray astronomy to operate on a satel- 
lite platform. A detailed description of the COMPTEL in- 
strument, which is sensitive to 7-rays in t he 0.8 -30 MeV 
range, can be found in Schonfelder et al. ( |1993| ). Briefly, 
the instrument consists of two planes of detector arrays, 
Dl and D2, separated by 1.58 m (see Fig. @). The Dl 
detector consists of seven cylindrical modules filled with 
NE 213A organic liquid scintillator. The D2 detector con- 
sists of 14 cylindrical Nal(Tl) crystals. The Dl scintillator 
material has a low average atomic number to optimize the 
occurrence of a single Compton scatter, while the D2 scin- 
tillator crystals have a high density and average atomic 
number to maximize their photon absorption properties. 
Each detector array is surrounded by a pair of overlapping 
anti-coincidence domes, manufactured of NE 110 plastic 
scintillator, to reject charged particle triggers of the tele- 
scope. The in-flight performance of the instrument is mon- 
itored with two calibration (CAL) units, each composed 
of a 60 Co-doped scintillator viewed by two 1/2 inch pho- 
tomultiplier tubes (PMTs), that provide tagg ed ph otons 
for in-flight energy calibration (Snelling et al. 
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COMPTEL instrument accepts and registers coincident 
triggers in a single D1-D2 module pair within the coinci- 
dence time window of ~ 40 ns in the absence of a veto sig- 
nal from the four charged particle shields as valid events. 
These interactions can be caused by a single photon or by 
multiple photons and/or particles. Among other param- 
eters, a time-of-flight (ToF) value and a so-called pulse 
shape discriminator (PSD) value in Dl are recorded for 
each event. The ToF is a measure of the time difference be- 
tween the triggers in the Dl and D2 detectors and is used 
to discriminate forward scattered (Dl— >D2) events, such 
as celestial photons with a ToF value of about 5 ns, from 
backward scattered (D2— »D1) background events which 
cluster around a ToF value of about — 5 ns. The PSD is 
a measure of the shape of the scintillation light pulse in 
the Dl detector. The energy loss characteristics of recoil 
electrons resulting from Compton scattering and of recoil 
protons resulting from neutron scattering are different, al- 
lowing one to reject many neutron induced events. The 
summed energy deposits in the two detectors, Ei + E 2 , 
are a measure of the total energy of the incident photon, 
Etot, while the photon scatter angle tp is determined from 
Ei and E2 through the Compton-scatter formula: 
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Fig. 1. A schematic view of the Com pton telescope 
COMPTEL (from Schonfelder et al. |l993j) . 



3. Characteristics of instrumental line 
background 

The particle and photon environment of COMPTEL pro- 
duces coincident interactions in the Dl and D2 detectors 
that pass all the logic and electronic criteria for a valid 
event. In the following, an overview is given of different 
types of instrumental background events in terms of their 
interaction process and location, as illustrated in Fig 



cos tp = 1 



Eo 



with m c 2 = 511 keV (1) 



This classification of background events (van Dijk pL996| 
provides a simple and versatile framework for discussing 
the COMPTEL instrumental background. ToF is of prime 
importance in identifying and rejecting instrumental back- 
ground events. A schematic representation of the ToF dis- 
tribution of valid events is depicted in Fig. |3[ 

• Event types A and B: events caused by the dou- 
ble scattering of a single photon. Any photon created in 
CGRO may produce this type of background event, which, 
if the scattering is from Dl to D2, is identical to a proper 
celestial event. It follows that the ToF distribution of for- 
ward scattered single photon events is identical to that of 
celestial photons and peaks around 5 ns in the ToF for- 
ward peak. Depending on where they originate, forward 
scattered single photon events may be rejected. In partic- 
ular, many of the photons that originate from below the 
Dl detector (such as type B events) can be eliminated by 
a selection on the scatter angle <p (see e.g. the event selec- 
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Fig. 2. An illustration of the main channels for the trig- 
gering of valid events due to instrumental background 
(adopted from van Dijk 1996 ). For simplicity, BATSE was 
omitted in this schematic view of CGRO. The various 
event types are explained in the text. 

tions in App. |a|). The remaining instrumental background 
from forward scattered single photons is therefore mostly 
due to photons originating in, around and above the Dl 
detector (such as type A events); type B events are negli- 
gible. Line emissions appear in energy space, particularly 
in E to t and E1-E2 space (see below). 

Single photons may also scatter in D2 before interact- 
ing in Dl. These backward scattered single photon events, 
however, are identifiable by their ToF distribution, which 
is confined in the backward peak at ToF values of about 
—5 ns. Finally, high-energy neutrons can undergo a double 
scattering process analogously to photons, however, many 
of these neutron events can be rejected by their PSD value. 

• Event types C and D: events caused by two or 
more photons both spatially and temporally correlated 
(so-called multiple photon or cascade events). In general, 
multiple photon events are more efficient in generating 
a background event than single photon events, since the 
probability for coincident interactions in both detectors 
increases with the number of emitted photons. Emission 
of two or more photons can occur from a small region 
on a time-scale shorter than the coincidence window of 
40 ns. Nucleons that have been excited above the first 
nuclear level, e.g. by proton or neutron interactions, may 
promptly emit a cascade of photons. The multitude of nu- 
clear excitation levels often results in a rather featureless 
continuum distribution in energy space. If, however, only 
a few transition levels are involved, characteristic features 
appear in E2 and E1-E2 space (see below). Multiple pho- 
tons may also arise from /3-particle bremsstrahlung, or the 
annihilation of a positron. Also, high-energy neutrons may 
induce the emission of one or more 7-ray photons outside 
of the Dl scintillator (thereby eluding rejection by PSD) 
and some may interact in the D2 detector. More compli- 
cated nuclear reactions such as the spallation (break-up) 
of a nucleus or the initiation of a shower of secondaries 
by an incident cosmic-ray particle or neutron may also 
produce multiple photons. 

Frequently, the emission of photons is nearly simulta- 
neous, i.e. on time-scales much shorter than the coinci- 
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Fig. 3. A schematic representation of the ToF distribution 
of valid events. Three major components can be discerned: 
the ToF backward peak and forward peak, centered at 
ToF values of about —5 ns and +5 ns, respectively, and 
an underlying continuum distribution. The backward peak 
is composed of all types of background events originat- 
ing in and around the D2 detector (these are not spec- 
ified individually in this illustration). The ToF forward 
peak contains the celestial signal as well as background 
events originating in and around the Dl detector, mostly 
of types A and C. The ToF continuum is dominated by 
background events of types D, E, and F originating in the 
instrument structure between the two detectors and the 
spacecraft structure in general. The relative magnitudes of 
the different components, which depend on E to t, are only 
represented approximately. 



dence window. The ToF value of these event types is de- 
termined by the location of the emitting nucleus relative to 
the Dl and D2 detectors. Multiple photons originating in 
the vicinity of the Dl detector (such as type C ev ents ) will 
peak slightly below a ToF value of 5 ns (see Sect. L2) and 
therefore contribute to the ToF forward peak. Photons 
emitted by nuclei in the spacecraft material between the 
Dl and D2 detectors (such as type D events) will inter- 
act in the two detectors near-simultaneously and have ToF 
values that are broadly distributed around zero, while pho- 
tons originating in and around the D2 detector will con- 
tribute to the ToF backward peak. The full, double-peaked 
ToF distribution of all multiple photon events anot only 
reflects the location of the emitting nucleus, but also the 
mass distribution of the entire spacecraft relative to the 
Dl and D2 detectors. 

• Event type E: events that are both spatially and tem- 
porally uncorrelated (the so-called random coincidences). 
The COMPTEL detectors are continuously exposed to a 
large flux of 7-ray photons. This inevitably leads to coinci- 
dent interactions that qualify as valid events. The photons 
producing these random coincidences are mostly of local 
or atmospheric origin. Since the two photons creating the 
event are unrelated to one another, and, in particular, not 
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correlated in time, these events are uniformly distributed 
in ToF. Type E events may also involve a neutron instead 
of a photon triggering the D2 detector, which has no PSD 
capability. 

• Event type F: events caused by two photons that are 
temporally correlated, but spatially uncorrelated. High- 
energy cosmic-ray particles or atmospheric neutrons may 
interact at several different locations along their path 
through CGRO. Individual interactions include those gen- 
erating type C and D events (e.g. spallations or showers). 
The whole interaction chain creating type F events is sim- 
ilar to multiple photon events, the main difference being 
that their ToF distribution depends on both the location 
of the interactions as well as the (relativistic) velocity of 
the primary particle. The ToF distribution of this type of 
event is broad and covers the entire coincidence window 
(also reflecting the spacecraft mass distribution). 

• Other processes, such as the interaction of a neutron 
in the D2 detector after producing a photon in the Dl 
system, or direct ionization losses of /3-particles created 
in the housings of the Dl scintillators, may also play a 
role. The event signatures of these and other, more com- 
plicated processes, however, will be similar to the event 
types described above. 

Event types A - D may arise from activation by pri- 
mary cosmic-ray protons or secondary particles as well as 
atmospheric neutrons. The time between the interaction 
of the proton or neutron within the telescope material 
and the actual triggering of a background event varies, 
since it depends on the decay time of the radio-isotopes 
produced. These processes and the resulting background 
events can be crudely separated into "prompt" and "de- 
layed" components. For prompt background events the 
time delay between the primary particle interaction in 
the instrument and the resulting emission of background 
photons is shorter than the coincidence window of 40 ns 
for the triggering of a valid event. Thus the intensity 
of prompt background components instantaneously fol- 
lows the (time- variable) incident local cosmic-ray flux (see 
Sect. [H]). For delayed events the time delay is longer than 
the typical length of the fast-logic veto signal of ~ 200 ns. 
In contrast to protons, neutrons can travel to any location 
in the spacecraft to produce 7-ray photons without trig- 
gering the veto system. Inside the veto domes, particularly 
in the Dl detector, prompt events can therefore only be 
produced by neutrons. Both protons and neutrons, how- 
ever, can produce delayed background events inside as well 
as outside the veto domes. 

As described in Sect. |[ events due to incident celestial 
photons have ToF values around 5 ns. Only background 
events with a similar ToF value will therefore interfere 
with astrophysical analyses. As illustrated in Fig. ||, a 
major portion of the instrumental background events, in- 
cluding those in the backward peak, can be eliminated by 
a ToF selection. The most important background event 
types in the ToF forward peak region around ToF= 5 ns 
are those originating in the Dl detector (types A and C), 



and some of the background events produced in the satel- 
lite structure (types D, E, and F). 

After ToF selection, the majority of the instrumental 
line background is expected to arise from activation of 
the Dl detector material because of the relatively high 
mass density and probability for triggering a background 
event as compared to the general spacecraft structure. The 
material composition of the Dl detector system therefore 
provides important clues as to which radioactive isotopes 
can effectively be produced and ultimately contribute to 
the instrumental line background in the low-Earth orbit 
of CGRO. The Dl support structure and the Dl mod- 
ule and PMT housings are mostly aluminium, the most 
abundant element in the instrument. The liquid scintil- 
lator NE 213A is composed of hydrogen and carbon. The 
quartz windows in the module housings contain silicon and 
oxygen. The PMTs and electronics boxes contain, among 
other elements, copper, nickel, and iron. 

Events from the instrumental line background produce 
conspicuous and characteristic features in energy space 
that can be exploited to distinguish them from the gen- 
eral (continuum) background. The ToF distributions of 
the line background and the continuum background are 
also different: the instrumental line background is concen- 
trated in the forward peak region, while the continuum 
background in energy is throughout the ToF continuum 
as well as the ToF forward peak. In the following, the 
characteristics of single photon (type A) events and mul- 
tiple photon or cascade (type C) events are illustrated 
with examples of (background) events from 2 D and 24 Na, 
respectively, both major contributors to the instrumental 
line background (see Sect. [I]). 

The instrumental 2.22 MeV photons emitted when 2 D 
is produced in the Dl scintillators (see Sect. 4.1) are sin- 
gle photon (type A) background events. The E to t, E2, Ei, 
and E1-E2 distributions of instrumental 2.22 MeV pho- 



tons for CDG event selections (described in App. A.l) as 
determined from Monte Carlo simulations are depicted in 
Fig. [|. The E to t distribution exhibits a peak at the energy 
of the primary photon, while the distributions in E2 and 
Ei are broad and relatively featureless. In E1-E2 space 
the event distribution of type A events follows the diag- 
onal Ei + E 2 = E 7 , with E 7 = 2.22 MeV for 2 D. The 
distributions in E to t and E1-E2 are the most important 
characteristics of single photon background events. 

The /3~-decay of 24 Na results in the emission of two 
photons with energies 1.37 MeV and 2.75 Me V causing 
type C line background events (see Sect. h2). The j3~- 
particle is of minor importance for the generation of a 
background event. Since it rarely escapes the support 
structure, it contributes to the background only through 
secondary bremsstrahlung photons. The 24 Na signature 
in E to t, E2, Ei, and E1-E2 as determined from Monte 
Carlo simulations is shown, for CDG event selections (see 
App. A.l), in Fig. [| The most important characteristic of 
24 Na, and any other isotope emitting multiple photons, is 
the E2 distribution, which exhibits peaks at the energies 
of each of the primary photons. In Ei the only conspicu- 
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Fig. 4. The E to t, E 2 , Ei, and Ei-E 2 distribution of instru- 
mental 2.22 MeV photons from 2 D production for CDG 
event selections as determined from Monte Carlo simula- 
tions. 
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Fig. 5. The E tot , E 2 , Ei, and Ei-E 2 distribution of instru- 
mental 24 Na events for CDG data selections as determined 
from Monte Carlo simulations. 



ous feature is the Compton edge of the 1.37 MeV photon 
(the Compton edge of the 2.75 MeV photon is suppressed 
by the CDG data selections). The E tot spectrum is more 
complex and exhibits less pronounced, line-like features 
just above the individual photon energies. These come 



from the absorption of one photon in D2 with the other 
photon scattering in Dl with an energy deposit near the 
Dl threshold. In addition, in E to t a shoulder is present 
at about 3.9 MeV. This is analogous to the sum peak in 
standard spectroscopy employing a single detector. Under 
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CDG data selections, this shoulder is due to the absorption 
of the 2.75 MeV photon in D2, with the 1.37 MeV photon 
scattering in Dl with an energy deposit at the Compton 
edge. The 24 Na event distribution in E1-E2 space is also 
very characteristic: the events cluster along two bands par- 
allel to the Ei axis located in E2 at the energies of the two 
decay photons. The band at 2.75 MeV in E2 extends in 
Ei from the threshold up to the Compton edge of the 
1.37 MeV photon. The band at 1.37 MeV in E 2 is sup- 
pressed by the CDG data selections. 

For other isotopes, if the energy of one of a pair of 
cascade photons is below the D2 threshold, then the two 
photons produce a cascade (type C) event only if the 
lower-energy photon scatters in Dl, and the higher-energy 
photon interacts in D2. The E2 and E1-E2 distribution of 
such a cascade photon pair is simple: there is only one 
photopeak in E2 and only one band in E1-E2. 

Radioactive decays that result in the (simultaneous) 
emission of two or more photons can produce line back- 
ground events other than of type C. For example, an indi- 
vidual photon of the emitted photon multiple can produce 
a type A event, provided all other photons escape from the 
instrument without interacting in any of the detectors. As 
apparent in the E1-E2 event distribution of 24 Na, however, 
the Etot signatures corresponding to the photon energies 
1.37 MeV and 2.75 MeV are much weaker than the cas- 
cade structure. Monte Carlo simulations show that less 
than 10% of the 24 Na background events are due to type 
A events. This illustrates the fact that multiple photon 
decays in the Dl detector material are more efficient in 
generating a background event than decays that result in 
the emission of a single photon only. 

4. Identified radioactive isotopes 

The major components of the COMPTEL instrumental 
line background can be attributed to eight individual iso- 
topes, namely 2 D, 22 Na, 24 Na, 28 Al, 40 K, 52 Mn, 57 Ni, and 



208 T1 (Weidenspointner 1999). Identification of these iso- 
topes was achieved in an iterative process, starting from 
the most prominent lines. The E to t and E2 spectra (see 
Fig. |^) are particularly useful for the identification of iso- 
topes that give rise to single photon (type A) or multi- 
ple photon (type C) events, respectively. The diagnostic 
power of Ei spectra is limited since they do not exhibit 
line features at the energy of the incident photons, but 
only rather broad features at the corresponding Compton 
edges. Due to the correlated signatures in E1-E2 space, in- 
dividual spectral features can be accentuated by applying 
suitable event selections. Viable isotope identifications are 
required to account self-consistently for spectral features 
in selected regions of the E1-E2 dataspace (see the detailed 
explanations in App. |^), as well as for their variation with 
time and/or incident cosmic-ray intensity (discussed in 
Sect. ||). The telescope response to individual isotopes was 
modelled through Monte Carlo simulations. The isotopes, 
their half-lifes, most important decay channels, and main 
production channels are summarized in Table |l| Below, a 
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Fig. 6. An illustration of the spectral distributions in Etot 
(top) and E2 (bottom) as a function of time. The two spec- 
tra plotted with thick lines represent the times of high- 
est (after the second reboost in May 1997) and lowest 
(from the beginning of the mission until the first reboost, 
May 1991 - Nov. 1993) contributions from long-lived back- 
ground isotopes. The three intermediate spectra, plotted 
with thin lines, cover the time periods Nov. 1993 - Oct. 

1994 (lowest), Oct. 1994 - Oct. 1995 (middle), and Oct. 

1995 - May 1997 (highest). 



review of the isotope identifications is given (a more de- 
tailed account on the COMPTEL instrumental line back- 



ground can be found in, e.g., Weidenspointner 1999) 



4.1. 2 D 

Thermal-neutron capture on hydrogen results in the pro- 
duction of a stable 2 D nucleus and the emission of a single 
2.22 MeV photon, as seen in E to t (see top panel in Fig. ||). 
The liquid scintillator NE 213A in the Dl detector mod- 
ules consists of 9.2% H and 90.8% C by mass, making 
it an efficient moderator for neutrons. The instrumental 
2.22 MeV photons are isotropically emitted from the Dl 
scintillators and constitute single photon (type A) back- 
ground events. The average moderation and absorption 
times in the Dl scintillator for an incident neutron are a 
few fis and about 3 x 10 -4 s, respectively, hence the instru- 
mental 2.22 MeV line is a quasi-prompt background com- 
ponent that follows the local, instantaneous cosmic-ray in- 
tensity (Weidenspointner et al. 1996). This quasi-prompt 
nature of the production/emission of the 2.22 MeV line 
is the reason for treating, in the following, the stable iso- 
tope 2 D as if it were unstable, and for referring to it as 
"short-lived". In principle, the initial scatterings of the 
incident neutrons could be identified by their PSD val- 
ues. However, the time-scale for the moderation and ab- 
sorption of an incident high-energy neutron is too long to 
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Isotope 


Half-Life 


Decay Modes and 
Photon Energies [MeV] 


Main 

Production Channels 


2 D 


prompt 


2.224 


1 H(ntbet,7) 


22 Na 


2.6 y 


P + (91%): 0.511, 1.275 
EC (9%) : 1.275 


27 Al(p,3p3n), 
Si(p,4pa;n) 


24 Na 


14.96 h 


/T: 1.37, 2.75 


27 Al(n,o), 
Al(p,3pn) 


28 Al 


2.2 min 


P~: 1.779 


27 Al(n th er,7) 


4o K 


1.28 x 10 9 y 


EC (10.7%): 1.461 


primordial 


52 Mn 


5.6 d 


EC (64%): 0.744, 0.935, 1.434 
13+ (27%): 0.511, 0.744, 0.935, 1.434 


Fe(p,x), Cr(p,x), 
Ni(p,x) 


57 Ni 


35.6 h 


/3+ (35%): 0.511, 1.377 

"C 1 ^ ^onO/\. 1 0*7*7 

-CjO (3U%): 1.377 


Ni(p,x), Cu(p,x) 


208 rpj 


1.4 x 10 10 y 
( 232 Th) 


P~ (50%): 0.583, 2.614 

p~ (25%): 0.511, 0.583, 2.614 


primordial 



Table 1. A summary of the iso- 
topes identified in the COMPTEL 
instrumental line background. For 
simplicity, only the photon energies 
of the most frequent decay modes 
are listed. If /3-decays are involved, 
the /3-particles have been included in 
the response simulations. The iden- 
tification of 208 Tl has to be consid- 
ered tentative. The label "prompt" 
for the half-life of the stable isotope 
2 D refers to the time-scale for the 
emission of the 2.22 MeV photon. 



associate the initial scatterings with the delayed neutron 
capture producing the 7-ray. The instrumental 2.22 MeV 
events therefore cannot be rejected by a selection on PSD. 

The bulk of the neutrons producing 2 D are expected 
to be of atmospheric origin, with secondary neutrons pro- 
duced in the spacecraft being of minor importance. This 
follows from measurements of the fast-neutron flux in 
the DI detector and Monte Carlo simulations of the pro- 
duction of secondary neutrons in cosmic-ray interactions 
(Morris et al. 1995a), as well as from an estimate for the 
neutron absorption efficiency of the Dl detector. 

4.2. 24 A/a 

The main production channels for 24 Na are neutron- 
capture reactions such as 27 Al(n,a) 24 Na and proton re- 
actions such as 27 Al(p,3pn) 24 Na. The Al structure of the 
Dl detector system is the primary source of 24 Na activa- 
tion. As described in Sect. spectral features of the 24 Na 
background in both the E to t and E2 spectra can be under- 
stood by considering its cascade nature (type C events). 
In particular, the two photons with energies of 1.37 MeV 
and 2.75 MeV produce line features at about 1.3 MeV and 
2.7 MeV in E 2 (see bottom panel in Fig. ||). 

In addition, the time-variation of the 24 Na event rate 
is consistent with expectations for an isotope with a half- 
life of 15 h that is produced during SAA passages (see 
Sect. 5.1). Also, ToF distributions summed using data se- 
lections that favour 24 Na events have the characteristics 
of type C background events, i.e. they peak slightly below 
a ToF value of 5 ns (see Fig. (?]) , corresponding to an aver- 
age distance of the location of the photon emission from 
the Dl module of about 20-30 cm. 

4.3. 40 K 

40 K is one of the contributors to the spectral line feature at 
about 1.4 MeV in E to t (see top panel of Fig. ||). A promi- 
nent line at this energy was present before launch in some 




ToF [nsec] 

Fig. 7. The measured ToF distribution of 24 Na events 
from the vicinity of the Dl detector. 



of the calibration data, and at that time was tentatively 
identified with primordial 40 K radioactivity, contained e.g. 
in the concrete of the buildings. The 1.4 MeV line was 
still present after launch, however, and did not vary with 
any orbital parameter, indicating an instrumental origin 
as well. 

A viable explanation for the instrumental 1.4 MeV line 
was t he potassium in the glass of the Dl PMTs (van Dijk 
1996 ). Electron capture by 40 K results in the emission of 
a single 1.46 MeV photon (type A event). This origin is 
supported by the fact that a disproportionate fraction of 
the photons at 1.46 MeV interact in the outer sections 
of the Dl scintillator, indicating that the photons enter 
the Dl detectors from the sides, consistent with the lo- 
cation of the Dl PMTs. Also, the observed ip distribu- 
tion at 1.46 MeV peaks at high (p values, again indicat- 
ing that the photons enter from the sides. These char- 
acteristics were reproduced in Monte Carlo simulations. 
Furthermore, the potassium content in the Dl PMTs de- 
rived from the observed 1.4 MeV line rate is consistent 
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with the manufacturer's specifications, as was confirmed 
by radiological measurements of a PMT in the laboratory. 
The mass of the PMT front window (made of Corning 
France 801.51) is about 24.9 g, with a potassium mass 
fraction of 6%; the mass of the PMT side and back glass 
(both made of Schott 8245) is about 40.7 g and 9.1 g, re- 
spectively, both with a potassium mass fraction of 0.14%. 
The front window contains more than 95% of the potas- 
sium in the PMTs and is closest to the detectors. The 40 K 
background therefore is dominated by the front window. 
The 40 K activity of the front window, normalized to its 
mass, is about 0.2 decays s _1 g -1 . 



4.4. 



22 



Na 



Routine data processing and the investigation of galac- 
tic 1.8 MeV line emission from radioactive 26 Al indicated 
the build-up of a broad background component at about 
1.5 MeV in E to t (see top panel of Fig. ||). Investigations of 
the energy distribution of this background feature, as well 
as of its ToF dependence, led to its identification with the 
/3 + -decay of 22 Na produced in the Dl detector material 



(Oberlack 1997 ). This background component is produced 
by a two photon cascade (type C event) with energies of 
~ 500 keV and ~ 1.3 MeV, consistent with the absorp- 
tion of the 1.275 MeV photon in D2, and the scattering 
of an annihilation 511 keV photon in Dl (the reverse pro- 
cess, involving absorption of the 511 keV photon in D2, is 
suppressed because 511 keV is below the D2 threshold). 
The ToF distribution of the 1.275 MeV 22 Na line in E 2 
is consistent with the ToF signature of the 24 Na cascade. 
The build-up of the feature is consistent with an isotope 
of half-life of 2.6 y (see Sect. |IJ ). 

The most important process for the production 
of 22 Na in the Dl detector are proton interactions 
in the aluminium structure, the main channel being 
27 Al(p,3p3n) 22 Na. Based on the production cross-sections, 
some 22 Na should also be produced in proton interactions 
with silicon, found in the glass of the Dl PMTs, such as 
Si(p,4pa:n). It is expected that most of the 22 N a produc- 
tion occurs during SAA passages (see Sect. 5.1). 



4.5. 28 AI 

The isotope 28 Al is produced by thermal-neutron cap- 
ture on aluminium, the most abundant element in the 
Dl detector. The existence of the strong 2 D line due to 
thermal-neutron capture on hydrogen indicates the pres- 
ence of a thermal-neutron flux in the Dl scintillators. The 
cross-section for thermal-neutron capture on aluminium 
has a value of 0.28 mb (c.f. 0.33 mb for hydrogen). Upon 
its /3~-decay, 28 Al can generate a type C background 
event: the 1.78 MeV photon is absorbed in D2, with a 
bremsstrahlung photon (from the f3~) scattering in Dl. A 
small fraction of the 28 Al decays result in type A back- 
ground events from the double scattering of the 1.78 MeV 
photon. 



From the beginning of the mission a weak, albeit sig- 
nificant, line feature at about 1.8 MeV has consistently 
been present in the E2 distribution (see bottom panel of 
Fig. ^). The position of this feature is independent of Ei, 
as expected for a type C event. Additional support for 
the identification of 28 Al comes from the fact that, as ex- 
pected, the 1.8 MeV feature is short-lived (see Sect. 5.2). 



4.6. 52 Mnand 57 Ni 

The above mentioned isotopes were not sufficient to ac- 
count entirely for the pronounced and broad ~ 1.4 MeV 
feature in E2 that built up over the mission. After sub- 
tracting the effects of the other isotopes, there was clear 
evidence for a build-up of the 1.4 MeV residual, whose 
half- life was estimated to be between that of 24 Na and that 
of 22 Na (assuming that it is due to a single isotope). The 
Ei dependence of the 1.4 MeV residual in E 2 suggested 
that the line is due to a /3 + -decay and therefore arises 
from type C background events. After detailed modelling 
of the E 2 and E to t spectra, taking into account the mate- 
rial composition of the Dl detector, it was concluded that 
the 1.4 MeV residual is due to a blend of two different 
isotopes, namely 52 Mn (1.43 MeV) and 57 Ni (1.38 MeV). 

The isotope 52 Mn is produced in SAA-proton interac- 
tions with the Fe, Cr, and Ni in the Dl detector^, found 
mainly in the electronics. Matters are complicated by the 
fact that 52 Mn can be produced in cither its ground state 
or an isomeric state. These two states have different half- 
lifes and decay schemes. The time variation of the 52 Mn 
event rate suggests that this isotope is more likely pro- 
duced in the ground state (Ti/2 = 5.6 d) than in the iso- 
meric state {T1/2 = 21.1 min), as described in Sect. 5.1. In 
this case, most 52 Mn events are due to the absorption of 
a 1.43 MeV photon in D2, accompanied by the scattering 
of another photon, from the radioactive decay or positron 
annihilation, in Dl. 

The isotope 57 Ni is expected to be mostly produced 
by proton interactions with Ni and Cu in the Dl detec- 
tor during SAA passages^. The 57 Ni background is mainly 
produced by the absorption of the 1.38 MeV photon in D2, 
with one of the 511 keV annihilation photons scattering 
in Dl (type C event). Similar to 52 Mn, the identification 
of 57 Ni, which has a half-life of about 36 h, is supported 
by modelling the observed long-term variation of its back- 
ground contribution (see Sect. 5.1). 



4.7. 



208 



Tl 



Similar to the 1.4 MeV feature, the 2.7 MeV feature in 
E 2 represents a blend of lines from more than a single 
isotope. In addition to 24 Na, at least one component with 
a very long half- life (exceeding that of 22 Na) is present. 



2 The Dl detector mass fraction of Fe, Cr, and Ni relative to 
Al is about 20%, 5%, and 6%, respectively. 

3 The Dl detector mass fraction of Cu relative to Al is about 
13%. 
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The position of the residual is about 2.6 MeV. 208 T1, which 
is part of the natural 232 Th decay chain, is the most viable 
candidate. 

The isotope 208 T1 undergoes /3~-decay through several 
channels, all of which involve the emission of at least two 
photons, implying that the 208 T1 multiple photon (type C) 
events are efficient in triggering the telescope. The half-life 
of 208 T1 is only 3.1 min. If the isotope is part of a natural 
decay chain, however, then its effective half-life is equal to 
the longest half-life of any of its parent isotopes, which in 
this case is the isotope 232 Th with a half-life of 1.4 x 10 10 y. 
The main channel for the production of background events 
by 208 T1 is the absorption of the 2.61 MeV photon in D2, 
with another photon scattering in Dl. 

It was assumed that 208 Tl is distributed like 40 K and is 
only present in the front glass of the Dl PMTs. This was 
tested by recording 7-ray spectra of a spare Dl PMT using 
a Ge spectrometer. Taking the COMPTEL results for the 
40 K and 208 T1 lines at face value, and assuming that both 
isotopes are equally distributed in the Dl PMTs, one ex- 
pects that the 208 T1 activity is about 10 -2 decays s" 1 g _1 
in the front window or about 6% of that of 40 K, based on 
the simulated efficiencies, corresponding to a 232 Th mass 
fraction of a few 10 -8 . Unfortunately, the laboratory mea- 
surement of 208 T1 was inconclusive. The 1.46 MeV 40 K 
line is detected at the 12a level above the general back- 
ground (1700 ± 140 counts). The strongest 208 T1 lines at 
2.61 MeV and 0.58 MeV, however, with intensities of 100% 
and 85.5%, were only detected at the 0.9a (50±56 counts) 
and 1.4a (128 ± 92 counts) level, respectively. Although 
consistent with expectations, these numbers are insuffi- 
cient to test the assumptions. Much longer integration 
times are needed. 

5. Variability of instrumental lines 

The activity of isotopes contributing to the COMPTEL 
background (except for the primordial radio-isotopes 40 K 
and 208 T1) in general varies in time. They are produced 
by interactions of primary and secondary cosmic-ray par- 
ticles and radiation-belt (SAA) particles within the in- 
strument material. The intensity of each of these particle 
species, and hence the level of activation, changes in time 
due to orbit details (e.g. the satellite's altitude, the geo- 
magnetic cut-off rigidity, the SAA radiation dose) and the 
solar cycle. The activity of a specific radioactive isotope 
then results from the competing processes of activation 
(production) and decay. 

The variation of the activity of a specific isotope is 
complex in general. If the isotope's half-life T±/ 2 and the 
typical time-scale(s) r for the intensity variation of the 
particle population(s) producing it are very different, how- 
ever, the activity variations can be described rather sim- 
ply. For T X j2 -C t, the isotope's activity approximatly 
follows the incident cosmic-ray particle intensity and its 
orbital variation (see Sect. 5.2). This is, e.g., the case for 
short-lived isotopes (T1/2 S; a few minutes) produced by 
primary cosmic-ray particles or neutrons, such as 2 D and 



28 



Al. For 7\/ 2 ^> t, a long-term build-up occurs since 
the activated nuclei do not decay away between consec- 
utive SAA transits (see Sect. 5.1). This is, e.g., the case 



for long-lived isotopes (T1/2 ^> 90 min, i.e. the orbital 
period) produced during SAA passages (which occur 6- 
8 times each day), such as 22 Na and 24 Na. Intermediate 
half-lifes produce more complicated variability. 

The long-term and orbital variations of spectral fea- 
tures are useful in identifying their physical origin. 
Spectral features arising from primordial radioactivity will 
show no variation. Prompt and short-lived components 
will vary with incident cosmic-ray intensity, i.e. on time- 
scales shorter than an orbital period. In addition, prompt 
and short-lived components may also vary over long time- 
scales if their production is sensitive to changes in the 
cosmic-ray intensity due to the solar cycle or the orbit 
altitude. Spectral features due to long-lived isotopes will 
exhibit variations on time-scales that reflect their half-life 
as well as the changes in SAA dosage. With regard to a 
spectral feature not yet identified, these long-term varia- 
tions provide a crude estimate of its half-life. With regard 
to a tentatively identified isotope, comparison of its time- 
dependent activity with an empirical model serves as an 
important cross-check for the correctness of its identifica- 
tion. 



5.1. Long-term variation 

The long-term variation of the activity of long-lived radio- 
nuclides arises from the combined effects of the isotopes' 
decay and the time history of the activation episodes dur- 
ing SAA transits. Activation outside the SAA by cosmic- 
ray particles is negligible (see e.g. Kurfess et al. 1989 ), 
since the average daily fluence due to SAA particles by 
far exceeds that of cosmic-ray particles. The count rate 
in the CAL-units can be used to model the production of 



long-lived radio- nuclides (Varendorff et al. 1997, hereafter 
referred to as "activity model"). The small plastic scin- 
tillators (thickness 3 mm, diameter 12 mm) of the CAL- 
units are sensitive to protons and other charged particles, 
but insensitive to secondary 7-ray photons. They there- 
fore provide a measure of the charged-particle flux at the 
instrument (Snelling et al. 1986 ). The trigger rates of the 
two CAL-units above both a low and a high threshold 
value are available for all times, including SAA transits. 
Most other instrument systems are switched off during 
SAA passages, including the veto shields which therefore 
cannot be used as charged particle monitors for the activ- 
ity model. The high-threshold trigger rate of CAL-unit B 
(tbh) is most useful because it is the least affected by noise 
and long-term efficiency degradation. Only telemetry data 
gaps interrupt this measure of the cosmic-ray intensity. 

To achieve continuous information on the charged- 
particle flux during every SAA passage, the activity model 
employs a neural net to describe tbh as a function of or- 
bit altitude, geographic longitude and latitude, time since 
launch (to include variations due to the solar cycle), and 
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orientation (azimuth and zenith) of the satellite relative 
to its velocity vector (to account for asymmetries in the 
incident SAA-proton flux, see e.g. Watts et al. 1989 ). The 
number of nuclei of a specific isotope as a function of time, 
N(t), is then given by: 



-\st 



+ N new (t) 



N(t + St) = N(t) 
and 

N new {t) = a ■ St ■ (jBYl(t) r BH,o) 



(2) 



where the time interval St usually is 16.384 s, A = In 2/2^/2 
is the decay constant of the isotope with half- life I\/2, a 
is the proportionality factor between the CAL-unit trigger 
rate and the isotope activation, and rBH,o is the value of 
tbh outside the SAA. The isotope activity A{t) then is 



A(t) = A • N(t) 



(3) 



The value of a is not well known as it depends on the 
cross-sections for the production of the isotope in the tele- 
scope material, particularly in the Dl detector. For each 
isotope, a scaling factor, which is the product of a and the 
efficiency for triggering a background event (see Sect. ^), 
must be determined from a measurement of the event rate 
from the COMPTEL data. 

A comparison of the measured event rates due to 
the long-lived isotopes 22 Na, 24 Na, 52 Mn, and 57 Ni with 
the predictions of the normalized activity model from 
May 1991 through July 1997 is depicted in the left panels 
of Fig. ^. In addition, the long-term variation determined 
for the event rates arising from the short-lived isotopes 
2 D and 28 Al, as well as from the primordial radio-nuclides 
40 K and 208 T1 are shown in the right panels of Fig. ||. 
These rates were obtained from the study of the galac- 
tic 26 Al 1.8 MeV line emission, which requires a mea- 
sure of the isotopes' background contributions for each 
observation period (see App. B.2) under imaging event 
selections (see App. A..2). The scatter of the eight sets of 
event rates, particularly those for the primordial radio- 
isotopes 40 K and 208 T1, is larger than what is expected 
from statistics alone. Part of the additional scatter is ex- 
pected from intrinsic inadequacies of the fitting procedure 
(such as the c orrec tion for the Earth-horizon angle selec- 
tion, see App. B^) and of the fit models employed (these 
cannot account for all structures in the E2 and E to t spec- 
tra, which particularly affects weak components, see below 
and App. |b|). Partly, however, the additional scatter may 
be attributed to known physical effects. Due to the pre- 
cession of the satellite orbit the radiation dose from SAA 
transits may vary considerably between observation peri- 
ods. In addition, the orientation of the satellite relative to 
the anisotropic SAA-particle flux affects the radiation dose 
in the Dl detector, which is the major source of instrumen- 
tal line background. The resulting differences in SAA radi- 
ation dose between observation periods contribute to the 
scatter in the event rates of the long-lived isotopes. Since 
the CAL-units are located half-way between the Dl and 
D2 detectors, shielding effects due to surrounding space- 
craft materials may result in less variation in tbh than in 



the Dl radiation dose, so that the activity model would 
underpredict this variation. Orientation effects may also 
result in variations of the effective cosmic-ray particle flux 
between observation periods and add some scatter to the 
event rates of short-lived isotopes. 

The long-term variation of the long-lived isotopes re- 
flects the variation of the SAA radiation dose, which de- 
pends, among other parameters, on orbit altitude and so- 
lar cycle. 24 Na is the best tracer of the changing inten- 
sity of the encountered SAA-particle fluxes, as it is the 
strongest component of the line background and has a 
half-live of about 15 h. From the beginning of the mission 
until the first reboost (around TJD 9280), the decay of 
the orbital altitude (see Fig. ||) results in a monotonic de- 
crease of the incident SAA-particle fluxes. The SAA radia- 
tion dose increased after the first reboost; the effect of the 
slow decrease of the orbital altitude is probably more than 
compensated by the solar cycle, which proceeded towards 
solar minimum, resulting in the observed net increase of 
SAA radiation dose over time. The second reboost (around 
TJD 10560) again resulted in a significant increase of ac- 
tivation during SAA passages. The long-term variation of 
the 52 Mn and 57 Ni event rates is similar to that of 24 Na, 
since the half-lifcs of these isotopes are of the same order. 
The long-term behaviour of the 22 Na event rate is quite 
different. Its 2.6 y half- life is significantly longer than the 
typical time-scale of the changes in SAA radiation dose. 
Therefore 22 Na never reaches an equilibrium between ac- 
tivation and decay, but is continuously built up, with the 
two reboosts resulting in increases in the slope of the 22 Na 
activity trend. 

The fact that the activity model reproduces the long- 
term variation of the event rates from long-lived isotopes 
confirms the model assumptions, in particular that ac- 
tivation predominantly occurs during SAA transits. The 
model can also verify the isotope identifications, as the 
activity of an isotope depends, among other parameters, 
on its half-life (see Eqs. || and ||). 

The background produced by the short-lived isotopes 
2 D and 28 Al is not expected to be influenced by SAA ra- 
diation dose, rather, any long-term trends will arise from 
variations in the average cosmic-ray intensity. The count 
rate in the 2.22 MeV line does not exhibit any significant 
long-term variation. The orbital variation of the event rate 
in the 2.22 MeV line (Weidenspointner et al. 1996| ) is sim- 
ilar to that of the fast-neutron flux (E„ > 12.8 MeV) 



in th e D l detector (Morris et al. 1995a ), as discussed in 
Sect. ^2. In addition, Morris et al. ( 1995a ) demonstrated 
that the solar cycle variation of the fast-neutron flux is 
much less than its orbit variation. It is therefore expected 
that the relatively weak solar cycle dependence of the 
2.22 MeV line is dominated by the much greater orbital 
variations, which presumably accounts for much of the ob- 
served scatter. Similar to 2 D, the isotope 28 Al is predomi- 
nantly produced by thermal-neutron captures. Hence 28 Al 
is expected to exhibit the same long-term variation as 2 D. 
However, the low count rate from 28 Al exhibits a slight, 
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Fig. 8. The measured event rates due to the long-lived isotopes 22 Na, 24 Na, 52 Mn, and 57 Ni, the short-lived isotopes 
2 D and 28 Al, and the primordial isotopes 40 K and 208 T1 as a function of time for the first 6 years (May 1991 through 
July 1997) of the CGRO mission for imaging event selections. In addition, the predictions of the normalized activity 
model for the long-lived isotopes is depicted with solid grey lines. 
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Fig. 9. The altitude of the CGRO orbit as a function 
of time from May 1991 through July 1997. The two re- 
boosts of the orbit around TJD 9280 (Oct.-Dec. 1993) and 
TJD 10560 (Apr.-May 1997) are clearly seen. The rapid 
decay of the orbit at the beginning of the mission results 
from the increased drag of the expanded atmosphere at 
solar maximum. 



increasing trend, hinting at deficiencies of the fit models 
and/or contributions from yet unidentified isotopes. 

The activity of the primordial radio-nuclides 40 K and 
208 T1 is constant over the duration of the mission. Both 
data sets are consistent with this fact, although a small 
long-term increase may be present for 208 T1. Again, this 



apparent increase may be due to deficiencies of the fitting 
procedure and/or due to yet unidentified line background 
components. 

The activity model can not only be used to predict the 
number of nuclei activated during SAA passages, but also 
to estimate the average daily SAA-proton fluence from the 
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daily increase in the number of proton-produced nuclei. 
For example, based on the measured 22 Na event rate a 
daily average SAA-proton fluence (E p > 100 MeV) o f 2.3 x 



10 protons cm was inferred by Varendorff et al. (1997) 0.02 



for the beginning of the mission during solar maximum at 
an altitude of 440 km. Considering the uncertainties in this 
measurement, and the large altitude gradient of the SAA- 
proton flux (Stassinopoulos 1989] ), this value is in good 
agreement with the prediction of 5 x 10 5 protons cm -2 for 
an altitude of 462 km (Dyer et al. 



1994) 



5.2. Variation with cosmic- ray intensity 

The prompt instrumental background closely tracks the 
local, instantaneous cosmic-ray intensity, which can, 
e.g., be parameterized by a geomagnetic cut-off rigidity. 
Another way of parameterizing the incident cosmic-ray 
intensity is to use the count rate of the anti-coincidence 
domes of the COMPTEL instrument, referred to as "veto 
rate" in the following^. To a good approximation, the 
prompt background components vary in proportion to the 
incident cosmic-ray intensity as monitored by the veto 
rate. The best example of this are 2 D and 28 Al ( Fig. |Io|) , 



determined for CDC event selections (see App. A.l ). In 
the following discussion, "veto rate" and "incident cosmic- 
ray intensity" are therefore used interchangeably. In con- 
trast to prompt background components, the activity of 
the primordial radio-isotopes 40 K and 208 T1 is indepen- 
dent of incident cosmic-ray intensity and hence does not 
vary with veto rate. The activity variation of long-lived 
isotopes such as 22 Na with veto rate is complex and de- 
pends on the isotopes' half-life as well as on the geophys- 
ical environment. Similar to the study of the long-term 
variation, the study of the variation of the event rate of a 
background component with cosmic-ray intensity can be 
helpful in identifying the responsible isotope. 

In the following, the physical origin of the complex 
and - in general - non-linear variation of the event rate 
of long-lived isotopes with incident cosmic-ray intensity 
as monitored by veto rate is illustrated with the isotope 
22 Na, which has a half- life of 2.6 y. To investigate the vari- 
ation of a line background component with veto rate over 
a given time period the events are sorted in energy spec- 
tra, one for each of a set of veto rate intervals, according to 
the veto rate value at the time the events were recorded. 
The event rate due to the background isotope under study 
is then determined for each veto rate bin in an iterative 



line fitting procedure described in App. B.l. Each of these 
rates represents the "average" event rate due to the iso- 
tope in a given veto rate interval over a given time period. 



4 This veto rate is the number of deadtime clock counts that 
occur during the presence of a veto signal. For each veto dome 
the deadtime clock counts are accumulated over the length 
of a telemetry packet (2.048 s) once during each superpacket 
(16.384 s) and recorded in the house-keeping datasets in units 
of [deadtime clock counts/2.048 s = DT] 
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Fig. 10. The veto rate variation of the count rate due the 
background isotopes 2 D and 28 Al as measured for CDG 
event selections. The dashed line is a linear fit to the mea- 
sured 2 D event rate. 



The veto rate variation of the event rate due to 22 Na 
over relatively short periods of time, during which the geo- 
physical environment did not change significantly, is illus- 
trated in Fig. [ll] for data from 1991-1993 (CGRO Phases I 
and II) and from 1994-1995 (CGRO Cycle 4). To a good 
approximation, the 22 Na event rate is independent of veto 
rate, as expected for a half-life of 2.6 years. The average 
22 Na event rate increased with time due to the build-up of 
this isotope during successive SAA passages (see Fig. |J). 
An increasing background contribution is also observed 
for long-lived isotopes with shorter half-lifes, such as 24 Na 
with Ti/2 = 15 h. However, due to a correlation between 
the encountered geomagnetic cut-off rigidity and the time 
since the last SAA transit, there is some correlation be- 
tween the event rate of these isotopes and veto rate, which 
is generally complex even over relatively short periods of 
time. Nevertheless, the following conclusions drawn from 
the veto rate variation of the 22 Na event rate apply to 
long-lived isotopes with shorter half-lifes such as 24 Na as 
well. 

The average variation of the event rate arising from 
22 Na over the extended time period from 1991-1996 
(CGRO Phase I through Cycle 5), during which the geo- 
physical environment changed significantly due to the first 
reboost (see Fig. ||) and the solar cycle, is depicted in 
Fig. [l^. Typically, a bump-like feature, hereafter referred 
to as the veto rate bump, appears at low veto rate values 
when studying the veto rate variation of the background 
contribution of long-lived isotopes over extended periods 
of time that include at least one reboost. This is even 
true for an isotope such as 22 Na, whose event rate is in- 
dependent of veto rate over short periods of time. It must 
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Fig. 11. The count rate from 22 Na, determined under 
CDG event selections, as a function of veto rate for data 
from 1991-1993 and from 1994-1995. 
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Fig. 12. The count rate from 22 Na, determined under 
CDG event selections, as a function of veto rate for data 
from 1991-1996. 

be emphasized, however, that changes in the geophysical 
environment do not always produce a bump-like feature, 
but rather result in a complex, non-linear variation of the 
contribution from long-lived background components with 
veto rate. 

To understand the origin of the veto rate bump it is 
necessary to realize that the changing geophysical environ- 
ment not only affects the activity of long-lived isotopes 
(see Fig. |J), but also the veto rate frequency distribu- 
tion, mostly due to variations of the incident cosmic-ray 
intensity with orbit altitude and solar cycle. This is il- 
lustrated in Fig. [l3[ depicting the normalized veto rate 
frequency distributions for the time periods 1991/1992, 
1994/1995, and 1997/1998 (i.e. CGRO Phase I, Cycle 4, 
and Cycle 7, respectively). For clarity, the distributions 
for other Phases/Cycles were omitted as they are fully 
consistent with the trend described in the following. The 
generic veto rate distribution is characterized by a pro- 
nounced peak at low veto rates, with a fiat tail extending 
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Fig. 13. The variation of the normalized veto rate 
frequency distribution with time, exemplified for veto 
dome 2. Depicted are the normalized veto rate frequency 
distributions, in units of 16.384 s time intervals (super- 
packets, SP) per veto rate bin, for the time periods 
1991/1992 (thin black line), 1994/1995 (grey line), and 
1997/1998 (thick black line). 
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Fig. 14. Illustration of the origin of the veto rate bump. 
Features such as the veto rate bump result from the com- 
bined effects of the time variation of the isotope activity 
(see Fig. 0) and of the veto rate frequency distribution 
(see Fig. O), as explained in detail in the text. 



to higher veto rate values. To a good approximation, the 
time variation of this generic veto rate frequency distribu- 
tion can be described by a slow but steady shift of the peak 
to higher veto rate values, with some slight changes of the 
shape of the distribution such as an increasing extension 
of the tail. 

As is illustrated in Fig. [l4|, the veto rate bump is noth- 
ing more than an artifact of the combined effects of the 
long-term variation in the activity of long-lived isotopes 
(see Fig. ||) and of the long-term shift of the veto rate 
frequency distribution (see Fig. For simplicity, let us 
consider a long-lived isotope whose activity is indepen- 
dent of veto rate over periods of time during which the 



14 



G. Weidenspointner et al.: The COMPTEL instrumental line background 



geophysical conditions do not change, such as 22 Na. Let 
us further consider a combination of two such time in- 
tervals of equal duration from early and late in the mis- 
sion. Late in the mission the isotope activity is higher 
than early in the mission, and the veto rate distribution 
is shifted to higher veto rate values. The two sets of data 
are represented by the two bars in Fig. |lj, with the grey- 
shading indicating the veto rate distribution in terms of 
the number of 16.384 s sampling intervals (superpackets, 
see Footnote 4) per veto rate interval, i.e. dark grey in- 
dicates that the corresponding veto rate values occurred 
more often (compare Fig. |l|). The average activity (deter- 
mined as described above) then increases from interval 1 
through interval 4 as the late-mission data become rel- 
atively more important. In interval 5 both sets of data 
contribute equally and the average activity is the mean 
of the individual activities. Finally, only data from the 
late mission contribute to veto rate interval 6. The net 
result for the average activity exhibits a general tendency 
to increase with increasing veto rate with a pronounced 
bump superimposed - similar to what has been measured 
for 22 Na (compare to Fig. |l2|). It can easily be demon- 
strated, following the above example, that the veto rate 
bump is an almost universal indicator for the presence of 
radio- isotopes with half-lifes ranging from several minutes 
to several years when combining data covering extended 
periods of time during which the geophysical environment 
inevitably changed. However, the veto rate bump cannot 
be used to estimate the half-life of an isotope, as can be 
done with the activity model. Also, it should be noted that 
non-linear variations of the event rate due to a long-lived 
isotope with veto rate other than the veto rate bump are 
to be expected once the isotope activity decreases and/or 
the veto rate frequency distribution shifts towards lower 
values as the mission continues. 

Among other things, a detailed understanding of the 
variation of the background contribution of long-lived iso- 
topes with veto rate is necessary for a reliable measure- 
ment of the cosmic diffuse gamma-ray background (CDG) 
at MeV energies. The COMPTEL analyses make use of 
the variable nature of the instrument background to iso- 
late the CDG signal, whic h is assumed to be constant 
(see e.g. Kappadath et al. 2000 , Weidenspointner et al. 
2001 ). In particular, above about 4 MeV the instrumental 
background is exclusively due to prompt processes. These 
are also present in the background at lower energies, but 
with large contributions from long-lived radioactive iso- 
topes. Because prompt background components vary, to a 
good approximation, linearly with incident cosmic-ray in- 
tensity as monitored by veto rate, they can be eliminated 
from the data by an extrapolation technique (Kappadath 
et al. 1996). In contrast, the event rate of long-lived iso- 
topes does not vary linearly with veto rate. The con- 
tributions of long-lived isotopes must therefore be sub- 
tracted, as a function of veto rate, before the veto rate 



lived background components, since only then can the 
event rate be assumed to vary linearly with veto rate. 
Subtraction of non-prompt background components after 
the veto rate extrapolation is only possible for constant 
background components, such as primordial radioactivity. 

6. Background line intensities: comparisons and 
implications 

The instrumental line background of an instrument is de- 
termined by its material composition, the radiation en- 
vironment, and the employed detection principle. In this 
respect it is interesting to compare the COMPTEL instru- 
mental lines with those identified in the data from other 
low-energy 7-ray experiments on satellite platforms such 
as the Nal Gamma-Ray Spectrometer (GRS) onboard the 
Solar Maximum Mission (SMM) and the Ge Gamma-Ray 
Spectrometer (GRS) flown on HEAO 3. In doing so, one 
must keep in mind that an instrument's energy resolu- 
tion is important for resolving lines. At 1 MeV, the en- 
ergy resolution (FWHM) of COMPTEL is 9.8% and 8.8% 
in E to t and E 2 , respectively (Schonfelder et al. |l993|) . At 
the same energy, th e res olution of GRS-SMM was about 
5.4% (Forrest et al. |l98l| ) and that of GRS-HEAO 3 was 
< 0.3% (Mahoney et al. |1980D . Therefore, the Ge spec- 
trometer onboard HEAO 3 is far superior for 7-ray line 
studies compared to COMPTEL and GRS-SMM, which 
have similar resolution. The CGRO satellite is kept in 
a circular orbit of 28.5° inclination at altitudes between 
about 330 km and 515 km (see Fig. |). The SMM and 
HEAO 3 missions were operated in circular orbits of 28.5° 
inclination at about 490-570 km altitude, and 43.6° incli- 
nation at about 500 km altitude, respectively. Similar to 
COMPTEL, the most abundant element in the material 
composition of both GRS-SMM (E. Chupp, priv. comm., 
1999) and GRS-HEAO 3 (Wheaton et al. |l98|) is Al in 
the instrument structures. Also, all instruments contain 
significant amounts of Fe, Ni, Cu, and Cr in, e.g., elec- 
tronics components. 

Five of the eight isotopes identified in the COMPTEL 
line background, namely 22 Na, 24 Na, 28 Al, 52 Mn, and 
57 Ni, are due to activation of Al and Fe, Ni, Cu, and Cr nu- 
clei in the instrument structure. As expected, each of these 
isoto pes has also been identified in GRS-S MM (S hare et al. 



1989[ ) and GRS-HEAO 3 (Wheaton et al. |989|) data. For 



extrapolation (Weidenspointner 199E). The validity of a 
veto rate extrapolation for removing prompt instrumental 
backgrounds thus depends on the absence of any long- 



GRS-SMM, however, an additional production channel for 
22 Na and 24 Na were proton and neutron reactions, respec- 
tively, on 23 Na in the Nal detectors. Other isotopes iden- 
tified in GRS-SMM and GRS-HEAO 3 that result from 
activation of the five elements and which have major lines 
below about 1 MeV should in principle also be present 
in the COMPTEL background, however, their events are 
suppressed by the Dl and D2 detector thresholds used 
in COMPTEL analyses. On the other hand, instrumen- 
tal lines due to activation of Cs and I (GRS-SMM) or Ge 
(GRS-HEAO 3) are absent in COMPTEL since they ex- 
hibit a ToF value not used in the analysis (I) or since these 
elements are not present in the instrument (Cs, Ge). 



G. Weidenspointner et al.: The COMPTEL instrumental line background 



15 



The different detection principles employed in the 
three instruments result in significant differences in the 
relative importance of specific isotopes, their detectabil- 
ity, and rejectability. COMPTEL requires coincident in- 
teractions in the Dl and D2 detector and therefore is 
particularly susceptible to multiple photon events such 
as the /3 + -decay of 22 Na or the /3~-decay of 24 Na (see 
Sect. ||). Consequently, multiple photon events are rela- 
tively more important than single photon events. On the 
other hand, individual photons from these decays produce 
photopeaks in E2 that can be used for identifying individ- 
ual isotopes and for determining their background con- 
tribution. In addition, the characteristic E1-E2 signature 
of multiple photon events provides a wide range of op- 
tions for their rejection (or enhancement). In the SMM 
and HEAO 3 spectrometers events are triggered by indi- 
vidual photons. Unlike the COMPTEL veto system, which 
consists of plastic scintillator domes, the anti-coincidence 
systems of the SMM and HEAO 3 spectrometers were 
made of Csl crystals and therefore sensitive to photons at 
7-ray line energies. For these two detectors therefore the 
probability for a radioactive decay to trigger the detectors 
is proportional to the number of emitted photons, as is the 
probability to trigger the anti-coincidence systems. To first 
order, the net probability of a radioactive decay for gener- 
ating a background event is therefore independent of the 
photon multiplicity for GRS-SMM and GRS-HEAO 3. As 
long as the decays occur outside the detectors, individual 
photons will give rise to photopeaks in the latter two in- 
struments, and these backgrounds are harder to suppress 
by event selections than in COMPTEL. As far as decays 
inside the detectors are concerned, /3-decays are partic- 
ularly hard to identify for GRS-SMM (e.g. activation of 
Na) and GRS-HEAO 3 (activation of Ge) since the added 
/3-continuum will broaden 7-ray lines beyond recognition. 
In COMPTEL, activation in the D2 detector is effectively 
eliminated by event (ToF) selection. 

Since COMPTEL is the first double-scattering 
Compton telescope operated in a near-Earth orbit, the 
instrumental background experienced during this mission 
may provide guidance for the design of future instruments. 
Below, some implications of the COMPTEL instrumental 
line background for conducting 7-ray line studies with this 
and future instruments are discussed. 

The by far strongest astrophysical 7-ray line signal de- 
tected by COMPTEL is 1.8 MeV line e missio n from 26 Al 
in the interstellar medium (Diehl et al. 1995a ). When ob- 
serving along the galactic plane, the average event rate 
due to this extended line source is about 7 x 1CT 4 s" 1 for 
imaging selections (see App. A. 2). Typical event rates due 
to point sources observed by COMPTEL in the light of 7- 
ray lines, such as the supernova remnant Cas A ( Ti at 
1.12 MeV, Iyudin et al. 1997), or the Type la supernova 
1991T ( 56 Co at 0.85 MeV and 1.24 MeV, Morris et al. 
1997bh , or the Vela region ( 26 A1 at 1.8 MeV, Diehl et al. 
1995b) are about 3 x 10 -5 s _1 . Typical event rates arising 



Isotope 


Efficiency 
(Imaging Sel.) 


Material 
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7 x 1CT 4 
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4 x 10~ 4 
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from individual background isotopes are about 10 s 
(see Fig. 0). Below 3 MeV, instrumental background lines 



Table 2. Background isotope properties relevant for as- 
trophysical 7-ray line studies. For each isotope the fol- 
lowing quantities are listed: the efficiency for producing a 
background event under imaging event selections, the ma- 
terials) in which the isotope is produced, and the mission- 
averaged activity. Further details are given in the text. 



acco unt for 10-50% of the total background rate (see 
Fig. p3.8| ). It follows that the signal-to-background ratio 
for astrophysical 7-ray lines in general is less than 1%, 
which can be enhanced to a few percent by fine-tuning 
the imaging event selections described in App. A. 2 for an- 
alyzing a specific 7-ray line. 

In this respect it is interesting to take a look at the effi- 
ciency, determined by Monte Carlo simulation, for various 
isotopes to produce a background event (see Table ||) . For 
imaging selections, this efficiency is about 7 x 10 -4 and 
2 x 10~ 4 for 2 D and 40 K, both of which give rise to type 
A (single photon) background events. A major source of 
type C multiple photon events is 24 Na with an efficiency of 
about 9 x 10 -4 . The efficiency for 208 T1, a minor source of 
type C events which is assumed to have the same spatial 
distribution as 40 K, is about 2 x 10 -3 . This is an order of 
magnitude larger than the efficiency for K and again il- 
lustrates that COMPTEL is more susceptible to multiple 
photon decays than to single photon decays. The more 
stringent CDG event selections reduce these efficiencies 
by factors of 4-7 and 2-3 for type A and type C back- 
ground events, respectively. For comparison, the detection 
efficiencies for celestial lines around 1-2 MeV are about 
2 x 10 -4 using imaging event selections. 

A goal for future Compton telescopes is to mini- 
mize background production by design. Passive material 
should be reduced, and manufactured from materials with 
low activation. Based on simulated efficiencies and the 
COMPTEL mass model, "mission-averaged" activities (in 
decays per g of material per s) were derived for the eight 
identified isotopes using data representing the first 7 years 
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of the mission (June 1991 to April 1998). These activities, 
together with the above mentioned simulated efficiencies 
and the activated materials, are summarized in Table 0. 
The relative yields for 52 Mn and 57 Ni were determined 
in hadron simulations (P. Jean, priv. comm., 1997). The 
numbers for the average activities are accurate within fac- 
tors of 2-3. Their major sources of uncertainty are the 
simulated efficiencies and the mass normalization, which 
are currently based on the assumption that the activation 
is homogeneous in a particular material. This may not 
be accurate for isotopes produced by thermal neutrons, 
in particular 28 Al. As a cross-check, the average activities 
were compared to simulations of activation of COMPTEL 
materials, neglecting the contributions of secondary parti- 
cles, for a solar-maximum radiation environment (P. Jean, 
priv. comm., 1997). The measured and the simulated ac- 
tivities are generally consistent within factors of a few for 
isotopes mainly produced by SAA protons. Hence these 
average activities may be used for estimating the line back- 
ground rates of future instruments in a similar orbit, tak- 
ing into account how the mass distribution differs from 
that of COMPTEL. 

7. Summary and discussion 

We have identified eight different isotopes in the 
COMPTEL instrumental line background, namely 2 D, 
22 Na, 24 Na, 28 Al, 40 K, 52 Mn, 57 Ni, and 208 T1. These iso- 
topes can account for the major instrumental background 
lines. Some minor instrumental lines, however, remain 
unidentified at this time. In addition, we have studied the 
variation of the event rate of these isotopes with time and 
incident cosmic-ray intensity, and determined the average 
activity of spacecraft materials. 

These results provide valuable insight into the physical 
processes that give rise to instrumental background lines 
in 7-ray detectors in low-Earth orbits. Because of their im- 
portance for MeV astronomy, these issues have repeatedly 
been studied (see, e.g., Dean et al. 1991 and references 
therein). Obviously, each detector has its own, unique in- 
strumental background. Nevertheless, background inves- 
tigations such as those reported in this work add to a 
growing pool of background experience to be found in the 
literature, which will prove to be important for the design 
of future instruments and the understanding of their data. 

The successful modelling of the time variation of the 
background contributions from the long-lived isotopes 
22 Na, 24 Na, 52 Mn, and 57 Ni in COMPTEL supports the 
conclusion that activation occurs predominantly during 
SAA passages. On average, CGRO passes through the 
SAA 6-8 times each day. In 1991 (at solar maximum and 
hence at minimum SAA-proton flux) the daily average of 
the incident flux of SAA protons with E p > 100 MeV 
was about 38 times higher than the corresponding inci- 



ample, during its mission the altitude of CGRO ranged 
from 350-500 km, which corresponds to a variation of the 
SAA-proton flux of a factor ~ 10, while the solar activity 
results in a variation by a factor ~ 2 at these altitudes 



(Stassinopoulos 1989). According to our activity model 
(Varendorff et al. 1997) the anisotropy of the SAA-proton 
flux at low satellite altitudes (Watts et al. 1989) has an 
important effect on the radiation dose and hence on acti- 
vation. These findings are consistent with earlier studies 
of GRS-SMM instrumental lines, in which activation due 
to (anisotropic) radiation-belt protons during SAA tran- 
sits has been identified as major source of (long-lived) line 
background and its variation with time (see Kurfess et al. 



1989, Share et al. 1989, and references therein). 



dent galactic cosmic-ray proton flux (Dyer et al. 1994), 
illustrating the importance of SAA passages with respect 
to the total deposited radiation dose. The SAA-proton 
flux varies strongly with altitude and solar cycle. For ex- 



In contrast, the activities of short-lived isotopes are 
mostly determined by the incident primary cosmic-ray and 
albedo-neutron fluxes. The enhanced production of these 
isotopes during SAA transits is of little importance for 
their background contribution outside the SAA, where as- 
trophysical data are recorded, because of their rapid decay, 
which precludes build-up. 

Two strategies for reducing instrumental line back- 
ground in future low-energy 7-ray detectors due to ac- 
tivation of long-lived isotopes present themselves. First, 
care should be taken to avoid detector materials which 
are easily activated by SAA protons (such as aluminium) 
or which carry primordial radioactivity. Second, instru- 
ments should be placed in orbits that minimize radiation- 
belt proton dosage, such as a low-altitude equatorial orbit 
below the radiation belts and outside the SAA, or a high- 
altitude orbit above the radiation belts (as scheduled for 
the INTEGRAL observatory). The flux of (low-energy) 
cosmic rays is significantly lower in a low-altitude orbit as 
compared to a high-altitude orbit - at the cost of a much 
higher albedo-neutron flux. Another strategy for limiting 
the impact of remaining instrumental lines in future in- 
struments is, e.g., to exploit the expected improvement 
in energy resolution to identify (and model) a large num- 
ber of instrumental background components. A compre- 
hensive treatment of the implications of the COMPTEL 
instrumental line (and continuum) background for future 
Compton telescopes is beyond the scope of this paper and 
will be given in a separate publication. 

Finally, although investigations of background lines are 
of interest in their own right, they are ultimately mo- 
tivated by the struggle to eliminate background events 
in astrophysical analyses. The COMPTEL instrumental 
background lines, particularly those from long-lived iso- 
topes, are a concern below about 4 MeV. At these en- 
ergies detailed understanding of the dataspace structure 
and variability of the instrumental lines is indispensable 
for optimizing many astrophysical investigations, such as 
that of the cosmic diffuse gamma-ray background (see e.g. 
Kappadath et al. 2000, Weidenspointner et al. 2001) or 
of the galactic 1.8 MeV line emission from 26 Al (see e.g. 
Oberlack 119971, Pliischke et al. ~ 



200C ) 



Comparing event rates due to 7-ray line sources and 
due to instrumental background it is evident that source 
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fluxes cannot be derived from global event rates with 
COMPTEL. This is only possible in imaging analysis, 
which exploits the characteristic cone-like distribution of 
source events in a three-dimensional data space spanned 
by the scatter angle (p and the direction (zenith and az- 
imuth) of the scattered photon (see Schonfelder et al. 
1993| ) . Instrumental background components exhibit com- 



paratively smooth distributions in this data space, allow- 
ing for their separation from the signal. Further exploita- 
tion of existing knowledge on the dataspace structure of 
instrumental lines will remedy some of the difficulties in 
the search for astrophysical 7-ray lines (see e.g. Morris et 
al. 1997b). Nevertheless, data selections that increase the 



signal-to-background ratio are clearly of great value. 

The importance of instrumental background suppres- 
sion/rejection can be illustrated with the original re- 
port of nuclear de-ex citation lines from the Orion region 
(Bloemen et al. 1994 ), which would have corresponded to 



an event rate of about 6 x 10 s spread over a rather 
broad energy range of 3-7 MeV. As was recognized in our 
recent analyses, the variations in the 24 Na event rate over 
the different observation periods of this galactic region 
are capable of producing spatial inhomogeneities in the 
background of that order, which then can be falsely at- 
tributed to a celestial signal. Only correlations of the sup- 
posed Orion signal with data space domains where 2 Na 
contamination is strongest provided the hints to suspect 
an instrumental background artifact as the origin of the 
claimed detection (Bloemen et al. 199E). 
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Appendix A: Event selections 

The results on the COMPTEL instrumental lines given in 
the main body of this paper have been derived using two 
different sets of event selections, described below. 



730 kcV < E 2 < 30 MeV 
800 keV < E tot < 30 MeV 

110 ch < ToF < 130 ch (A.l) 
PSD low (Er) < PSD < PSD high (E 1 ) 
6° < (p < 38° 

C + <fi = £ < 40° 

where (p is the photon scatter angle (see Eq. |l|) , and £ is the 
angle between the direction of the scattered photon and 
the telescope z-axis. The selection on £ defines a field-of- 
view of half opening angle 40°. The values of PSDi ow and 
PSD h i gh are 40 ch and 110 ch at 70 keV, 55 ch and 95 ch 
at 1 MeV, 55 ch and 92 ch at 10 MeV, and 40 ch and 90 ch 
at 20 MeV. At intermediate Ei energies, these boundaries 
are linearly interpolated. In addition, the horizon of the 
Earth, as seen by the instrument, is required to be at least 
5° outside the ficld-of-view. 



A. 2. Imaging analysis 

In COMPTEL imaging analyses such as that of the 26 Al 
1.8 MeV line emission from the galaxy, 7-ray albedo pho- 
tons from the Earth's atmosphere are rejected by deter- 
mining for each event the minimal angular distance be- 
tween its possible origins in the sky and the Earth hori- 
zon. This angular distance, denoted e below, is required 
to exceed a minimum value (see e.g. Schonfelder et al. 



1993 ). Hence, the effective held-of-view in imaging analy- 



sis changes with time as the instrument orbits the Earth. 

The standard set of imaging event selections used in 
this paper are as follows: 



70 keV < Ei < 20 MeV 
650 keV < E 2 < 30 MeV 
750 keV < E tot < 30 MeV 

115 ch < ToF < 130 ch (A.2) 
ch < PSD < 110 ch 
0° < (p < 50° 
e > 5° 



A.l. CDG analysis 

In the COMPTEL analysis of the cosmic diffuse gamma- 
ray background (CDG), the intense atmospheric 7-ray 
background is eliminated by selecting data from times 
when the Earth was sufficiently far outside a circular held- 
of-view defined by event selections. Detailed descri ptions 
of these selections can be fou nd in, e.g., Kappadath ( 1998 ) 
and Weidenspointner ( 199S ) . 

The standard set of CDG event selections used in this 
paper are as follows: 

70 keV < Ei < 20 MeV 



Appendix B: Line fitting procedures in analyses 
of CDG and 26 AI 

Many astrophysical analyses, such as that of the cos- 
mic diffuse gamma-ray background (CDG), of the galactic 
1.8 MeV line emission from 26 Al, or of 7-ray line emission 
from supernovae, can be optimized by detailed modelling 
of the instrumental line background. Below, the (similar) 
procedures used to determine the background event rates 
due to individual isotopes employed by Weidenspointner 
et al. ( |200lD , and Oberlack ( |1997[ ) and Pliischke et al. 
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( |2000D in investigations of the CDG| and of the 26 Al emis- 
sion, respectively, are described. 

B.l. CDG analysis 

Accounting for the event rates due to the instrumental 
line background, in particular its long-lived components, 
is important for measuring the CDG at MeV energies (e.g. 
Kappadath et al. 2000| , Weidenspointner et al. 2001). As 
explained above, the background contributions of long- 
lived isotopes have to be subtracted before prompt and 
short-lived backgrounds can be eliminated by veto rate 
extrapolation. The event rates due to the eight identi- 
fied background isotopes producing the major instrumen- 
tal lines are determined, as a function of veto rate, in an 
iterative procedure by fitting a set of three E2 spectra and 
one E to t spectrum for each veto rate interval]^ 

The rationale of the iterative fitting procedure is illus- 
trated in Figs. [B.l| and |Q . The Ei-E 2 distributions of 
the eight isotopes as obtained from Monte Carlo simula- 
tions are schematically depicted for CDG event selections 
in Fig. |B.l| . In general, there is considerable overlap in the 
E1-E2 distributions of individual isotopes (in particular 
around 1.3 MeV in E2), which precludes an independent 
determination of the isotopes' background contributions. 
Therefore, an iterative procedure was introduced, which 
starts at the highest energies in Ei and E2, where am- 
biguities are minimal, and then proceeds down to the in- 
creasingly complex structures at lower Ei and E2 energies. 
The Ei-E 2 ranges of the three E 2 spectra and the E to t 
spectrum, chosen such as to enhance or suppress individ- 



ual lines or spectral features, are indicated in Fig. B.2. 
The Ei-E 2 ranges covered by the second and third E 2 fit 
and by the E to t fit overlap, hence the results of these fits 
are not statistically independent. The overlap is caused by 
the Etot fit, which was introduced to properly separate the 
background events from 2 D and 28 Al. Inclusion or omis- 
sion of the Etot fit therefore is a trade-off between sys- 
tematic and statistical uncertainty. By iteratively fitting 
the second and third E2 spectrum and the E to t spectrum, 
both the systematic and the statistical uncertainty in the 
2 D and 28 Al event rates are minimized (other isotopes are 
hardly affected by the overlap, see below). In additon, the 
iterative approach ensures the self-consistency of the de- 
termined isotope background contributions. Also included 
in Fig. B.2 arc lines of constant (p and E to t- Comparison of 



5 A similar procedure consisting of four differ ent fit s was uti- 
lized in the CDG analysis of Kappadath et al. (200C) to deter- 



mine the background contributions from the five isotopes ( 2 D, 
22 Na, 24 Na, 28 Al, and 40 K (fixed)) identified at that time. 

6 Ideally, the event rates due to these isotopes would be de- 
termined by fitting the two-dimensional event distribution in 
E1-E2. However, at present such two-dimensional fits are not 
feasible for two reasons. First, in general the statistics is not 
sufficient in the CDG analysis, and second, we do not yet have 
a reliable model of the E1-E2 distribution of the continuum 
background, which has a complex shape as it is a blend of 
various components. 
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Fig. B.3. An example of a fit of the first E2 spectrum, 
which is used to determine the event rate due to 24 Na. In 
addition to the total fit, the 24 Na template (solid line), the 
two unidentified features (dashed and dash-dotted lines) , 
and the exponential continuum are indicated. 



Figs. 



B.l 



and 



B.2 provides a first indication of the many 
options for fitting individual lines, which can be enhanced 
or suppressed through the choice of the fit regions, and in 
addition through selections on E to t and/or (p. In particu- 
lar, event selections may be used to suppress unidentified, 
long-lived spectral features, which cannot be eliminated 
by veto rate extrapolation (see Kappadath et al 
Weidenspointner et al. 2001). 



2000 



The E1-E2 ranges of the three E2 spectra and the E to t 
spectrum indicated in Fig. B.2 were chosen for the itera- 
tive fitting procedure for the following reasons. The first 
E2 spectrum, covering the E1-E2 region of 950-1250 keV in 
Ex and 2000-3500 keV in E 2 (see Figs. [BT] and [B]|), allows 
us to determine the event rate from 24 Na. The signal from 
this isotope is optimized by selecting Ei energies around 
the Compton edge of the 1.37 MeV photon interacting in 
Dl and around the photopeak of the 2.75 MeV photon 
interacting in D2 (see Fig. |). The 2 D event rate is deter- 
mined from fitting the Etot spectrum (Etot 1810-2 800 k eV, 
Ei 70-950 keV, and E 2 730-2800 keV, see Figs. SO and 



B.2|) . The second E 2 spectrum (Ei 500-950 keV, E 2 1500- 
3500 keV, see Figs. B.l and B.2) is used to determine the 
event rate from 28 Al. Finally, the third E 2 spectrum is in- 
tended for determining the background contributions from 
the /3+-decays of 22 Na, 52 Mn, and 57 Ni, with the 270- 
350 keV range in Ei being optimized for the Compton 
edge of 511 keV photons, and the E2 range covering the 
energies 1100-3500 keV (see Figs. |y] and |bJ). To opti- 
mize the signal of the instrumental lines, which originate 
in the Dl detector material, a ToF range of 2.5-7.5 ns was 
selected for the spectra (compare Fig. ||) . 

The three E2 spectra and the Etot spectrum are an- 
alyzed in an iterative procedure consisting of eight fits, 
listed below. The contributions from the primordial radio- 
nuclides 40 K and 208 T1 are not determined from the fits, 
but calculated from their known ( 40 K, see Sect. 4.3) or 
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Fig. B.l. A schematic representa- 
tion of the simulated E1-E2 dis- 
tributions of the eight isotopes 
identified in the instrumental line 
background for CDG event se- 
lections. The "vertical" and "di- 
agonal" bands represent multi- 
ple photon and single photon 
events, respectively (compare to 
the detailed individual diagrams in 
Sect. 0). 



Fig. B.2. An illustration of the 
E1-E2 ranges of the three E2 spec- 
tra and the Etot spectrum used 
to determine the event rates due 
to the background isotopes. Also 
plotted are lines of constant ip 
(solid) and E to t (dashed). 



estimated ( 208 T1, see Sect. 4.7) activities based on Monte 
Carlo simulations. Similarly, once the background contri- 
bution of an isotope has been determined from any spec- 
trum, the isotope's contribution to any other spectrum 
can be predicted based on Monte Carlo simulations. The 
eight identified isotopes can account for the major instru- 
mental lines, however, some weak lines or spectral fea- 
tures remain unidentified at this time. In the fits, some 
of these unidentified lines were described by gaussians to 
minimize systematic errors in the determination of the 
event rates of the identified components. The unidentified 
components are genuinely different from those identified 
since their variation with cosmic-ray intensity, as well as 
their dependence on event parameter selections, are dif- 



ferent (see Kappadath et al. 2000, Weidenspointner et 
al. 200 ip . The identified isotopes are represented by tem- 
plates obtained from Monte Carlo simulation. These tem- 



plates have not been smoothed as smoothing inevitably 
increases the systematic uncertainty due to distortions of 
the tem plate sha pe. The small "spikes" in the templates 
in Figs. B.3-B.6 are an artifact of the plotting software. 
The eight fit steps are: 

1. E2 Fit 1: The event rate due to the isotope 24 Na is 
determined (see Fig. B.3 ). 

2. E to t Fit: The 2 D event rate is estimated; the contribu- 
tions from 24 Na and 208 Tl are fixed. All other isotopes 
are neglected. 

3. E 2 Fit 2: The 2 D and 28 Al event rates are estimated; 
the contributions from 24 Na and 208 T1 are fixed. 

4. E2 Fit 3: The background contributions due to the 
/3 + -decays of 22 Na, 52 Mn, and 57 Ni are estimated; the 
contributions from 2 D, 24 Na, 28 Al, 40 K, and 208 T1 are 
fixed. 
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Fig. B.4. An example of a fit of the E to t spectrum for 
determining the count rate in the 2 D 2.22 MeV line. In 
additon to the total fit, the 2 D template (solid line), an 
unidcntifcd component at about 2.3 MeV (dashed line), 
the continuum contributions from 24 Na and 208 T1 (solid 
and dash-dotted lines), and the exponential continuum 
(dashed line) are indicated. 

5. Etot Fit: The 2 D event rate is determined; the contri- 
butions from all other isotopes ( 22 Na, 24 N a, 28 Al, 40 K, 
52 Mn, 57 Ni, and 208 T1) are fixed (see Fig. 

6. E2 Fit 2: The event rate due to the isotope 28 Al is de- 
termined; the contributions from 2 D, 24 Na, and 208 T1 



are fixed (see Fig. B.5). 

7. E 2 Fit 3: The event rates due to the /3 + -decays of 22 Na, 
52 Mn, and 57 Ni are determined; all other isotopes are 
fixed (see Fig. |B.6| ) . 

8. E tot Fit: The 2 D event rate is re-determined as in Fit 5. 

In this procedure, Fits 5-7 re-iterate Fits 2-4 to ensure 
convergence, which is tested by comparing the results of 
Fits 5 and 8. In general, the two fits yield nearly identi- 
cal results. This self-consistent determination of the event 
rates of long-lived isotopes inevitably requires the deter- 
mination of all isotope contributions, whether long-lived, 
short-lived, or prompt. In the CDG analysis, subtraction 
of the contributions from 2 D, 24 Na, and 28 Al are based 
on Fits 8, 1, and 6, respectively (see Weidenspointer et 
al. 2001). The event rates due to the /3 + -decays of 22 Na, 
52 Mn, and 57 Ni are derived from the results of Fit 7. As 
pointed out before, the contributions from the primordial 
radio-nuclides 40 K and 208 T1 need not be determined from 
the line fits, but can be computed based on the known or 
estimated activities. 
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Fig. B.5. An example for a fit of the second E2 spectrum, 
which is used to determine the event rate from 28 Al. In 
addition to the total fit, the templates for 2 D (the line 
feature at ~ 1.6 MeV, fixed), 28 Al (the 1.78 MeV line), 
24 Na (the strong line at 2.75 MeV, fixed) and 208 T1 (the 
weaker line at 2.61 MeV, fixed), the three unidentified 
features (dashed, dashed-dotted, and dashed lines), and 
the exponential continuum are indicated. 
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Fig. B.6. An example for a fit of the third E2 spectrum, 
which is used to determine the event rates from 22 Na (dot- 
ted line), 52 Mn (thick solid line), and 57 Ni (thin solid line). 
Also depicted are the fixed components ( 2 D: solid line at 
1.9 MeV, 24 Na: dash-dotted line, 28 Al: weak solid compo- 
nent below 1.9 MeV, 40 K: dashed line, 208 T1: solid line at 
2.6 MeV), as well as the total fit, the exponential contin- 
uum, and the unidentified 2.93 MeV feature. 



B.2. 26 Al analysis 

The analysis of the galactic 1.8 MeV line emission from 
26 Al is, similar to that of the CDG, affected by the instru- 
mental line background and its temporal variation. The 
galactic 1.8 MeV line emission is determined in the en- 
ergy band 1.7-1.9 MeV, using adjacent energy bands for 
constructing a model for the background in this so-called 



line interval (comprehensive descriptions of this approach 
can be found in, e.g., Knodlseder 1997 and Oberlack 1997). 
In particular, the scatter angle {0) distribution of the 1.7- 
1.9 MeV background model is derived from an interpola- 
tion of the (p distributions in narrow adjacent energy in- 
tervals (1.6-1.7 MeV and 1.9-2.0 MeV). Due to the long- 
term variation of the instrumental line background (see 
e.g. Fig. @), the ratio of the number of counts in the line 



G. Weidenspointner et al.: The COMPTEL instrumental line background 



21 



1.02 

■B LOOP 
co 
tr 

c 0.98 hTi 

" 0.96 : 



0.92 t 





i 








3H 




H 


Mi 


. rf-Uti I . 1: 



8500 



9000 



9500 10000 
Time [TJD] 



10500 



11000 



Fig. B.7. Count ratios (number of counts in the line inter- 
val divided by the number of counts in the adjacent energy 
intervals) as determined for observations at galactic lat- 
itude \b\ > 40°. The ratios, and a linear model of their 
time variation, are shown before (grey) and after (black) 
subtraction of the instrumental line background. 



interval and in the adjacent energy intervals is decreas- 
ing with time (see Fig. B.7), mostly due to the build-up 
of 22 Na, a major component in the 1.6—1.7 MeV band in 
Etot (Oberlack 1997). To eliminate the time dependent 



contamination of the background reference from adjacent 
energies, the number of counts due to each component of 
the instrumental line background has to be determined for 
each individual observation period, employing a procedure 
outlined below. After subtraction of the instrumental line 
background, the ratio of counts in the line interval and 
the adjacent energy inter vals to a good approximation is 
constant in time (see Fig. B.7), and data from individual 
observation periods can be summed to analyze the galactic 
1.8 MeV line emission. 

The procedure used in the 26 Al analysis to determine 
the background contributions of the eight identified back- 
ground isotopes is a modified version of the CDG proce- 
dure described in App. 



B.l 



The modifications are moti- 
vated by differences in the event selections applied in the 
two analyses, particularly differences in the elimination 



of atmospheric background (see Apps. A.l and A. 2). One 
of the consequences is that in the 26 Al imaging analysis 
events with larger Cp values are accepted than in the CDG 
analysis. For some isotopes, notably 24 Na, this results in 
significant changes of the E1-E2 distribution as compared 
to that for CDG event selections. Consequently, the op- 
timal fit regions in E1-E2 space are somewhat different. 
Also, in contrast to the CDG analysis, in the 26 Al anal- 
ysis the (p distribution of the accepted events is different 
for each observation period due to the orbit dependent re- 
jection of atmospheric background. The simulated energy 
distributions for each background isotope therefore have 
to be corrected for the specific (p distribution of each indi- 
vidual observation period (Oberlack 1997). The correction 
for selections to reject atmospheric background has been 
calculated assuming a homogeneous illumination of the Dl 
detector. Since there is a certain edge enhancement in the 
illumination of the Dl modules for background produced 
in the Dl structure, slight differences in these corrections 
result in small distortions of the templates, and thus add 
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Fig. B.8. An example for an E to t fit in the 26 Al analysis. 
The upper panel depicts the Etot spectrum before (line 
histogram) and after (data points) subtracting the contri- 
butions from the identified background isotopes (indicated 
by smooth curves with different line types). The middle 
panel shows the E to t spectrum after isotope subtraction, 
together with the continuum background fit. The bottom 
panel gives the residuum of the fit. 



to the observed scatter in the determined isotope event 
rates (see Fig. ^). As each observation period is indepen- 
dently analyzed for its instrumental line background con- 
tamination, the background modelling in the 26 Al analysis 
yields as a welcome byproduct the long-term variation of 
the event rates due to the eight identified background iso- 
topes depicted in Fig. ||. 

The four steps of the iterative line fitting procedure in 
the 26 Al analysis are similar to those in the CDG anal- 
ysis, hence we content ourselves with describing only the 
differences. In the first step, two E1-E2 regions are fitted 
simultaneously to obtain the 24 Na event rate (Oberlack 
1997): the Ei-E 2 region depicted in Fig. |B.2|, and an ad- 



ditional E1-E2 region extending from 2000 to 2700 keV 
in Ei and from 1100 to 2000 keV in E 2 . Analogous to 
the first E1-E2 region, the additional E1-E2 region op- 
timizes the 2 Na signal from simultaneous interactions of 



22 



G. Weidenspointner et al.: The COMPTEL instrumental line background 



Fit 1 



it 20- 



Fit 2 



2 3 

r. 






60 




50 




40 


u 








0) 


30 


cr 
0) 


lH 






20 




10 









Fit 3 



20 - 
15 - 
10 - 

5 - 



Fit 4 



2 3 



2 3 

r. 



Fig. B.9. The distribution of the xt values for each of the 
four fit steps in determining the instrumental line contam- 
ination in the 26 Al analysis (Pliischke et al. 2000|) . 



the 1.37 MeV and 2.75 MeV photons in D2 and Dl, respec- 
tively, which are no longer suppressed due to the larger ac- 
cepted (p values for imaging selections. The E1-E2 region 
for the second step extends down to 700 keV in Ei, and is 
used for determining the 28 Al rate and for obtaining start 
parameters to determine the event rates of the remain- 
ing isotopes, particularly 22 Na, in Fit 3. However, unlike 
in the CDG analysis, the contributions of the primordial 
isotopes 40 K and 208 T1 are not fixed in the 26 Al analy- 
sis, but determined in Fit 3 and 1, respectively. As far as 
determining the contributions of each individual isotope 
is concerned, the procedure ends here, since it has been 
demonstrated in the CDG analysis that one fit cycle is 
sufficient to obtain self-consistent isotope rates. This ad- 
vanced background treatment provided the basis for the 
latest COMPT EL 26 Al all-sky maps (e.g. Oberlack |l~997 , 
Pliischke et al. p000[ ). 

To investigate the extent to which the identified iso- 
topes can account for spectral features in Etot, an ad- 
ditional fit is performed (see Fig. |B.g| ). In this E to t fit 
(Ei > 70 keV, E 2 > 650 keV), the contributions from 
the eight background isotopes are fixed at the values ob- 
tained from the three previous E2 fits; only the continuum 
background, modelled by a power law with an exponen- 
tial turn-over at low energies, is varied. As can be seen, 
the major spectral features are accounted for by the eight 
identified isotopes, however, some minor features remain 
unidentified at this time. 

The distribution of the \ 2 V values for each of the four fit 
steps, as obtained in the analysis of Pliischke et al. ( 2000D , 
is shown in Fig. B.9. The quality of Fits 1 and 3, in which 
the event rates of most of the background isotopes are 



determined, is acceptable. The xt distribtion for the E to t 
fit, however, indicates that further improvement of the 
instrumental line background modelling is still possible. 
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